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ABSTRACT 

The formation of protostellar disks out of molecular cloud cores is still not fully understood. Under 
ideal MHD conditions, the removal of angular momentum from the disk progenitor by the typically 
embedded magnetic field may prevent the formation of a rotationally supported disk during the 
main protostellar accretion phase of low mass stars. This has been known as the magnetic braking 
problem and the most investigated mechanism to alleviate this problem and help removing the excess 
of magnetic flux during the star formation process, the so called ambipolar diffusion (AD), has been 
shown to be not sufficient to weaken the magnetic braking at least at this stage of the disk formation. 
In this work, motivated by recent progress in the understanding of magnetic reconnection in turbulent 
environments, we appeal to the diffusion of magnetic field mediated by magnetic reconnection as an 
alternative mechanism for removing magnetic fiux. We investigate numerically this mechanism during 
the later phases of the protostellar disk formation and show its high efficiency. By means of fully 3D 
MHD simulations, we show that the diffusivity arising from turbulent magnetic reconnection is able 
to transport magnetic flux to the outskirts of the disk progenitor at time scales compatible with the 
collapse, allowing the formation of a rotationally supported disk around the protostar of dimensions 
~ 100 AU, with a nearly Keplerian profile in the early accretion phase. Since MHD turbulence is 
expected to be present in protostellar disks, this is a natural mechanism for removing magnetic flux 
excess and allowing the formation of these disks. This mechanism dismiss the necessity of postulating 
a hypothetical increase of the Ohmic resistivity as discussed in the literature. Together with our 
earlier work which showed that magnetic flux removal from molecular cloud cores is very efficient, this 
work calls for reconsidering the relative role of AD for the processes of star and planet formation. 
Subject headings: diffusion — ISM: magnetic fields — MHD — turbulence — star formation — 
accretion disks 



1. INTRODUCTION 

Circumstellar disks (with typical masses ~ 0.1 Mq 
and diameters ~ 100 AU) are known to play a funda- 
mental role in the late stages of star formation and also 
in planet formation. However, the mechanism that al- 
lows their formation and the decoupling from the sur- 
rounding molecular cl oud core progenitor is still not fully 
understood (see, e.g.. iKrasnopolskv et al.l[2011l for a re- 
cent comprehensive review). Former studies have shown 
that the observed embedded magnetic fields in molecu- 
lar cloud cores, which imply magnetic mass-to-flux ra- 
tios re lative to the critical value a few tinies lar ger than 
unity (|Crutcheil2005HTroland fc Crutcheil2"00l are high 
enough to inhibit the formation of rationally supported 
disks during the main protostellar accretion phase of low 
mass stars, provided that ideal MHD applies. This has 
been known as the magnetic braking problem (see e.g. , 
Allen et al.l[200l IGalli et aL| [200l [Price fc Batd [20071: 



Hennebelle fc Fromang|l2008riMellon fc Lill2008D . 

Proposed mechanisms to alleviate this problem and 
help removing the excess of magnetic flux during the 
star formation process include non-ideal MHD effects 
such as ambipolar diffusion (AD) and, to a smaller 
degree, Ohmic dissipation effects. The AD, which 
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was fi rst discussed in this context by [Mestel fc Spitzeii 
([19561). has be en extensively investig ated si nce then (e.g., 
i Spitzed [19681 [Nakano fc Tademanil [19721: [Mo uschoviasl 

. '1979: Naka no fc Nakamura[[1978r Shu 198a 

3r fc Mouschoviasl '1992, 1993 



Lizano fc Shu' 1989; Fiedler 



Li et a l. 2008; Fatuzzo fc Adams 2002; Zwcibcl 2002:). 



In principle, AD allows magnetic flux to be redis- 
tributed during the collapse in low ionization regions as 
the result of the differential motion between the ion- 
ized and the neutral gas. However, for realistic lev- 
els of core magnetization and ionization, recent work 
has shown that AD does not seem to be sufficient to 
weaken the magnetic braking in order to allow rota- 
tionally supported disks to form. In some cases, the 
magnet ic br aking has been found to be even enhanced 
bv AD dMel lon & Li 2009; Krasnooolskv & Konigl 2001 



Basu fc Mouschovias 1995; Hosking fc Whitw orth i2004[ : 
Duffin fc Pudritz 2009; Li et al. 2011). ^ These findings 
motiv ated iKrasnopolskv et al.i ([2010) (see also [Li et al.l 
[201 It ) to examine whether Ohmic dissipation could be ef- 
fective in weakening the magnetic braking. They claimed 
that in order to enable the formation of persistent, ro- 
tationally supported disks during the protostellar mass 
accretion phase a highly enhanced resistivity, or "hyper- 
resistivity" rj > 10^^ cm^s~^ of unspecified origin would 

See however a recent work that investigates the effects of AD 
in the triggering of magneto-rotational instability in more evolved 
cold, proto- planetary disks wh ere the fraction of neutral gas is 
much larger llBai £~Stonell2011l) . 
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be required. Although this value is somewhat depen- 
dent on the degree of core magnetization, it implies 
that the required resistivity is a few orders of magnitude 
larg er than the classic microsc opic Ohmic resistivity val- 
ues (|Krasnopolskv et aLll2010f). 

On the other hand . iMachida et afl (120101) (see also 
llnutsuka etahl 120101: IMachida et al.l 120 111) performed 
core collapse three-dimensional simulations and found 
that, even with just the classical Ohmic resistivity, a 
tinny rotationally supported disk can form at the be- 
ginning of the protostellar accretion phase (see also 
IDapp fc Basull20100 and grow to larger, 100- AU scales 
at later times. They claim that the later growth of the 
circumstellar disk is caused by the depletion of the in- 
falling envelope. As long as this envelope remains more 
massive than the circumstellar disk, the magnetic brak- 
ing is effective, but when the circumstellar disk becomes 
more massive, then the envelope cannot brake the disk 
anymore. In their simulations, they assume an initially 
much denser core than in Krasnopolsky et al. work, 
which helps the early formation of a tiny rotating disk 
facilitated by the Ohmic diffusion in the central regions. 
But they have to wait for over 10"" yr in order to allow a 
large-scale rotationally supported, massive disk to form 
(see discussion in §5). While this question on the ef- 
fectiveness of the Ohmic diffusion in the early accretion 
phases of disk formation deserves further careful testing, 
we here discuss an alternative more efficient mechanism 
to diffuse the magnetic flux based on turbulent reconnec- 
tion. 

Before addressing this new mechanism, it is 
crucial to note first that the concep t of "hyper- 
resistivity" previously me nt ioned (s ee alsolStrausslll986l: 
Bhattachariee fc Hameiril 119861 : iDiamond fc Malkovl 
2003f ) is not phvsicallv ju stified and therefore one cannot 
rely o n it (see criticism in lKowal et al.l[2"009HEvink et al.l 
120111) . Therefore, the dramatic increase of resistivity is 
not justified. 

However, one may notice that the problem we face with 
the magnetic field is not of dissipation of magnetic flux, 
where, indeed, ordinary resistivity is necessary, but of 
magnetic diffusion. Frequently in the literature, mag- 
netic diffusion is disregarded for highly conductive astro- 
physical flows. This is based on the concept of magnetic 
field being frozen-in into conductive fluids and thus mov- 
ing together with the fluid. However, the frozen-in con- 
dition is true in the absence of magnetic reconnection. 
When the latter is present, changes in the magnetic field 
topology are allowed; 

On the basis of iLazarian fc Vishniad (I1999D (hence- 
forth LV99) model of reconnection. iLaza rianl (|2005D 
(henceforth L05) suggested a way of removing mag- 
netic flux in the process of star formation, present- 
ing the concept of "reconnection d iffusion" of magnetic 
field in turbulent media (see also ILazarian fc Vishniad 
|2009() . This concept w as appl ied successfully to star for- 
m ation in [S antos-Lima, etaLl (120 10) ( sec also discussion 
in lLazarian et al.ll2010tlde Gouveia D al Pino et al. 2011, 
and n ew numerical calculation s in Leao et al. 2011). 
There, iSantos-Lima et al.l (|2010[ ) (henceforth SXIO) per- 
formed 3D MHD simulations of turbulent flows and 
found a decoupling between the gas and the magnetic 
field due to reconnection under one-fluid approximation, 
i.e., without having ambipolar diffusion. In addition, in 



the presence of gravity, they found a decrease of the mag- 
netic flux-to-mass ratio with increasing gas density in the 
center of the gravitational potential well, both for sys- 
tems starting with equilibrium distributions of gas and 
magnetic field and for dynamically unstable, collapsing 
systems. This provides evidence that the process of mag- 
netic flux removal by turbulent reconnection diffusion is 
important to quasi-static subcritical clouds and also to 
collapsing supercritical cores. Thus it is natural to ex- 
plore the consequences of the process in other systems. 

In this work, we investigate this effect on the removal 
of ma g,netic flux from collapsing, rotating protostellar 
cores. ILazarian fc Vishniad ()2009D argued that the re- 
moval of magnetic fle lds from c i rcums tellar disks re- 
ported in the work of iShu et al.l (|2006D is due to pro- 
cesses of turbulent reconnection, but have not provided 
any quantitative study of the effect. A preliminary 
discussion of the quantita ti ve re sults can be found in 
Ide Gouveia Dal Pino et all ()2011|) . Here, we show by 
means of 3D MHD numerical simulations that turbu- 
lent magnetic reconnection diffusivity enables the trans- 
port of magnetic flux to the outskirts of the core at time 
scales compatible with the collapse time scale, thus allow- 
ing the formation of a rotationally supported protostclar 
disk with nearly Keplerian proflle. 

In §2, we discuss the theoretical foundations of our 
work, in §3 we describe the numerical setup and ini- 
tial conditions, in §4 we show the results of our three- 
dimensional (3D) MHD turbulent numerical simulations 
of disk formation, and in §5, we discuss our results within 
a bigger picture of reconnection diffusion processes. Our 
summary is presented in §6. 

2. THEORETICAL CONSIDERATIONS 

The magnetic diffusion mechanism that we address is 
the process deeply rooted in microscopic physics of how 
magnetic fields behave in highly conductive flows. The 
textbook way to characterize these flows is to use the 
Lundquist number S = LxVa/ti, where rj is Ohmic dif- 
fusivity and Lx and Va are the scale of system at hand 
and the Alfven velocity, respectively. For astrophysical 
systems is very large and therefore S" 3> 1. The brute 
force numerical study of such systems is not feasible as 
the corresponding Lundquist numbers of the numerical 
simulations are much smaller. However, if the Lundquist 
number of the flow does not control magnetic reconnec- 
tion this opens prospects of modelling magnetic diffusion 
in astrophysical systems. 

The classical magnetic reconnection that follows the 
textbook Sweet-Parker scenario (see Figure 1, upper 
panel), depends on the Lundquist number and it is slow 
for astrophysical systems. Indeed, in the model all the 
matter moving with the speed Vrec over the scale 
should be ejected with the Alfven velocity through a thin 
slot A. The disparity of astrophysical typical scales 
and the scale A determined by microphysics, i.e. resis- 
tivity, makes Sweet-Parker reconnection rate negligibly 
small for most of astrophysical applications, including 
the case of accretion disks. However, this scenario is not 
valid in the presence of ubiquitous astrophysical turbu- 
lence (see Figure 1, lower panels). For the turbulence 
case, LV99 showed that the reconnection becomes inde- 
pendent of the resistivity, i.e. becomes fast, as the out- 
flow region A gets determined by magnetic fleld wander- 
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Fig. 1. — Upper plot: Sweet— Parker model of reconnection. The 
outflow is limited by a thin slot A determined by Ohmic difl^usivity. 
The other scale is an astrophysical scale Lx 2> A. Middle plot: 
reconnection of weakly stochastic magnetic field according to LV99. 
The outflow is limited by the magnetic field wandering. Low plot: 
an in dividual small-scale reconnection region. From iLazarian et al.l 



ing. This challenges the well-rooted concept of magnetic 
fie ld frozenness for t he case of turbulent fluids (see more 
in lEvink et al.lf20Tll henceforth ELVll) and provides an 
interesting way of removing magnetic flux out of, e.g. 
accretion disks . The LV99 model was successfully tested 
numerically in lKowal et al.l ()2009D . 

The justification of the reconnection diffusion con- 
cept in LOS i s based on the LV99 model (see also 
ILazarian et al.l 2004 for the case of magnetic reconnec- 
tion in partially ionized gas). Numerical effects are al- 
ways a concern when dealing with reconnection and mag- 
ne tic field diffu s ion. Indeed, unlike the numerical tests 
m iKowal et al.l (|2009l) . in our simulations the reconnec- 
tion events are happening on small scales where numer- 
ical effects are important. Precisely because of that, 
the numerical expe riments with anomalous resistivity in 
IKowal et al.l ()2009( ) are of key importance. There, using 
a numerical setu p with high resoluti on in a magnetic re- 
connection layer, IKowal et al.l ()2009[ ) showed that in the 
presence of turbulence the local non-linear enhancements 
of resistivity were not important. This confirmed the cor- 
responding analytical prediction in LV99 (see more dis- 
cussion in ELVll). Appealing to that finding, we claim 
that the reconnection diffusion that we observe in our 
simulations is a real effect and not a numerical artefact*. 
Analytical studies summarized in ELVll also support the 
notion that magnetic fields are generically not frozen-in 

As it was discussed in LV99 and shown even more explicitly 
in ELVll, the plasma effects that can enhance the local reconnec- 
tion speed are not important in the presence of turbulence which 
induces magnetic field wandering. 



when conductive fluids are turbulent. In view of them, 
the concept of reconnection diffusion in LOS looks very 
natural. 

As emphasized before, the concept above of magnetic 
fiux transport in turbulent flows by reconnection dif- 
fus ion has been already suc cessfully tested numerically 
bv lSantos-Lima et al.l (|2010f ) for idealized models of star 
forming clouds. In the present paper, we study whether 
the same concept can entail substantial changes for the 
magnetic fleld removal in the formation of protostellar 
disks. An extension to accretion disks in general can be 
also foreseen (see below). 

3. NUMERICAL SETUP AND INITIAL DISK CONDITIONS 

To investigate the formation of a rotationally sup- 
ported disk due to turbulent reconnection magnetic fiux 
transport, we have integrated numerically the following 
system of MHD equations: 

|^ + V.(pu) = (1) 

p + '^^ " = -CsV/3+^(V X B) X B-pVP + f 

(2) 

dA 

— = uxVxA - 77oftmV X V X A (3) 

where p is the density, u is the velocity, ^I* is the 
gravitational potential generated by the protostar, B is 
the magnetic field, and A is the vector potential with 
B = V X A H- Bea;t (where ^ext is the initial uniform 
magnetic field), f is a random force term responsible 
for the injection of turbulence. An isothermal equation 
of state is assumed with uniform sound speed Cs- 

We solved the MHD equations above in a three- 
dimensional domain using a second-order shock- 
capturing Godunov scheme and second-order Runge- 
Kutta time integration. We employ ed a modified version 
of the code originally developed bv IKowal et al.l ([2007), 
using the HLL Riemann solver to obtain the numerical 
fiuxes in each time step. 

In order to compare our results with those of 
iKrasnopolskv et al.l ()2010[ ) , we have considered the same 
initial conditions as in their setup. 

Our code works with cartesian coordinates and vector 
field components. We started the system with a collaps- 
ing cloud progenitor with initial constant rotation (see 
below) and uniform magnetic fleld in the z direction. 

Given the cylindrical symmetry of the problem, we 
adopted circular boundary conditions. Eight rows of 
ghost cells were put outside a inscribed circle in the xy 
plane. For the four outer rows of ghost cells, we adopted 
fixed boundary conditions in the radial direction in every 
time step, while for the four inner ghost cells, linear in- 
terpolation between the initial conditions and the values 
in the interior bound of the domain were applied for the 
density, velocity, and vector potential. With this imple- 
mentation the vector potential A has kept its initial null 
value. Although this produces some spurious noisy com- 
ponents of B in the azimuthal and vertical directions in 
the boundaries, these are too far from the central regions 
of the domain to affect the disk evolution. In the z direc- 
tion, we have applied the usual open boundary conditions 
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(i.e., zero derivatives for all conservative quantities: den- 
sity, momentum and potential vector). We found this 
implementation far more stable than using open bound- 
ary conditions in the x and y directions, or even in the 
radial direction. Besides, adopting circular rather than 
square boundaries prevented the formation of artificial 
spiral arms and corners in the disk. 

For modelling the accretion in the central zone, the 
technique of sink particles was i mplemented in the cod e 
in the same way as described in iFederrath et al.l (|2010l ) . 
A central sink with accretion radius encompassing 4 cells 
was introduced in the domain. The gravitational force 
inside this zone has a smoothing spline func tion identical 
to that presented in IFederrath et all ()201Cl[ ). We do not 
allow the creation of sink particles elsewhere, since we 
are not calculating the self-gravity of the gas. We note 
that this accreting zone essentially provides a pseudo in- 
ner boundary for the system and for this reason the dy- 
namical equations are not directly solved there where ac- 
cretion occurs, although we assure momentum and mass 
conservation. 

The physical length scales of the computational do- 
main are 6000 AU in the x and y directions and 4000 
AU in the z direction. A sink particle of mass 0.5 Mq 
is put in the center of the domain. At t = 0, the gas 
has a uniform density po — 1.4 x 10""'^^ g cm"'^ and a 
sound speed of Cg — 2.0 x lO"* cm s^^ (which implies 
a temperature T w 4.8/1 K, where Jt is mean molecu- 
lar weight in atomic units). The initial rotation profil e 
IS = Cs tanh(i?/i?c) (as in iKrasnopolskv et al.ll2010l ). 
where R is the radial distance to the central z-axis, and 
the characteristic distance Rc = 200 AU. 

We employed a uniform resolution of 384x384x256 
which for the chosen set of parameters implies that each 
cell has a physical size of 15.6 AU in each direction. The 
sink zone has an accretion radius of 62.5 AU. ^ 

Although we are interested in the disk that forms in- 
side a radius of approximately 400 AU around the central 
axis, we have carried out the simulations in a much larger 
region of 6000 AU in order to keep the dynamically im- 
portant central regions of the domain free from any outer 
boundary effects. 

4. RESULTS 

We performed simulations for four models which are 
listed in Table 1. Model hydro is a purely hydrody- 
namical rotating system. All the other models have 
the same initial (vertical) magnetic field with intensity 
Bz = 35 /xG. In order to have a benchmark, in the model 
named resistive we included an anomalous high resis- 
tivity, with a magnitude about 3 orders of magnitude 
larger than the Ohmic resistivity estimated for the sys- 
tem, i.e., r) = 1.2 x 10^° cm ^ s~^. According to the 
results of IKrasnopolskv et al.l ()2010D . this is nearly the 
ideal value that the magnetic resistivity should have in 

^ We note that in the two-dimensional simulations of 
IKrasnopolskv et al.l 12010) . they use a non- uniform mesh with a 
maximum resolution of 0.2 AU in the central region. The employ- 
ment of a non-uniform mesh could be advantageous in this prob- 
lem, allowing a better resolution close the protostar. However, in 
the present work since we are dealing with turbulence injection in 
the evolving system, the use of a uniform mesh has the advan- 
tage of making the effects of numerical dissipation more uniform 
and therefore, the analysis of the turbulent evolution and behavior 
more straightforward. 



order to remove the magnetic flux excess of a typical col- 
lapsing protostar disk progenitor and allow the formation 
of a rotationally sustained disk. We have thus included 
this anomalous resistive model in order to compare with 
more realistic MHD models that do not appeal to this 
resistivity excess. 

We have also considered an MHD model with turbu- 
lence injection (labeled as turbulent model in Table 1). 
In this case, we introduced in the cloud progenitor a 
solenoidal turbulent velocity field with a characteristic 
scale of 1600 AU and a Mach number Ms w 4 — 5 increas- 
ing approximately linearly from t — Q until t = 3 x IQ^'^ 
s (or « 3000 yr). These parameters result an estimated 
turbulent diffusivity which is of the order of the anoma- 
lous diffusivity employed in resistive model: rjturb ^ 
VturbLinj ~ 10^° cm^ s^^. The induced turbulent ve- 
locity field has been intentionally smoothed beyond a ra- 
dius of 800 AU, by a factor exp{-[i?(A[/) - 800]V4002}, 
in order to prevent disruption of the cloud at large radii. 
The injection of turbulence was stopped at t = 4.5 x 10^-^ 
s ~ 0.015 Myr. From this time on, it naturally decayed 
with time, as one should expect to happen in a real sys- 
tem when the physical agent that injects turbulent in the 
cloud ceases to occur. 

The last of the models (which is labeled as ideal MHD) 
has no explicit resistivity or turbulence injected so that in 
this case the disk evolves under an ideal MHD condition. 



TABLE 1 

Summary of the models 



Model 


Bo (mG) 


VOhm 

(cm^ s"-"-) 


Vturb 
(cm^ S~^) 


hydro 











resistive 


35 


1.2 X 10^0 





turbulent 


35 





~ 1020 


ideal MHD 


35 









Figure 2 shows face-on and edge-on density maps of 
the central slice of the disk for the four models at « 0.03 
Myr. The arrows in the top panels represent the direction 
of the velocity field, while those in the bottom panels 
represent the direction of the magnetic field. 

The pure hydrodynamical model in the left panels of 
Figure 2 clearly shows the formation of a high density 
torus structure within a radius « 300 AU which is typical 
of a Keplerian supported disk (see also Figure 3, top-right 
panel) . 

In the case of the ideal-MHD model (second row panels 
in Figure 2), the disk core is much smaller and a thin, 
low density outer part extends to the outskirts of the 
computational domain. The radial velocity component 
is much larger than in the pure hydrodynamical model. 
The bending of the disk in the core region is due to the 
action of the magnetic torques. As the poloidal field lines 
are dragged to this region by the collapsing fluid, large 
magnetic forces develop and act on the rotating flow. 
The resulting torque removes angular momentum from 
the inner disk and destroys its rotational support (see 
also Figure 3, upper panels). 

The third column (from left) of panels in Figure 2 
shows the resistive MHD model. As in Model 1, a torus 
(of radius ~ 250 AU) with a rotationally dominant ve- 
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locity field is formed and is surrounded by a flat, low 
density disk up to a radius of ~ 500 AU. Compared to 
the ideal MHD model (second column), the structure of 
the magnetic field is much simpler and exhibits the fa- 
miliar hourglass geometry. 

The last column (on the right of Figure 2) shows the 
ideal MHD model with injected turbulence (labeled tur- 
bulent). A high density disk arises in the central region 
within a radius of 150 AU surrounded by turbulent de- 
bris. From the simple visual inspection of the velocity 
field inside the disk one cannot say if it is rotationally 
supported. On the other hand, the distorted structure of 
the magnetic field in this region, which is rather distinct 
from the helical structure of the ideal MHD model, is 
an indication that magnetic fiux is being removed by the 
turbulence in this case. The examination of the velocity 
and magnetic field intensity profiles in Figure 3 are more 
elucidative, as described below. 

Figure 3 shows radial profiles of: (i) the radial velocity 
vn (top left), (ii)the rotational velocity (top right); 
(iii) the inner disk mass (bottom left); and (iv) the ver- 
tical magnetic field Bz (bottom right) for the models of 
Figure 2. vr and w$ were averaged inside cylinders cen- 
tered in the protostar with height h — 400 AU and thick- 
ness dr = 20 AU. Only cells with a density larger than 
100 times the initial density of the cloud (po — 1.4x 10~^^ 
g cm~^) were taken into account in the average evalua- 
tion. The internal disk mass was calculated in a similar 
way, but instead of averaging, we simply summed the 
masses of the cells in the inner region. The magnetic 
field profiles were also obtained from average values in- 
side equatorial rings centered in the protostar with radial 
thickness dr = 20 AU. 

For an ideal rotationally supported disk, the centrifu- 
gal barrier prevents the gas to fall into the center. In 
this ideal scenario, the radial velocity should be null (at 
distances above the accretion sink zone). The top left 
panel of Figure 3 depicts the curves of the radial veloci- 
ties for the four models. Above the central sink accretion 
zone {R > 62.5 AU), the hydrodynamical model (hy- 
dro) is the prototype of a rotationally supported disk, 
the radial velocity being smaller than the sound speed 
(cs = 2 X 10"* cm s~^) inside the formed disk. In the 
ideal MHD model, the effect of the magnetic flux braking 
partially destroys the centrifugal barrier and the radial 
(infall) velocity becomes very large, about three times 
the sound speed. The MHD model with anomalous re- 
sistivity instead, shows a very similar behavior to the 
hydrodynamical model due to the efficient removal of 
magnetic flux from the central regions. In the case of 
the turbulent model, although it shows a persisting non- 
null radial (infall) speed even above R > 62.5 AU this is 
much smaller than in the ideal MHD model, of the order 
of the sound speed where the disk forms (between the R 
« 70 AU and « 150 AU). 

The top right panel of Figure 3 compares the rotational 
velocities w$ of the four models with the Keplerian pro- 
file vk = \/ GA'h/R. All models show similar trends to 
the Keplerian curve (beyond the accreting zone), except 
the ideal MHD model. In this case, the strong suppres- 
sion of the rotational velocity due to removal of angular 
momentum by the magnetic field to outside of the inner 
disk region is clearly seen, revealing a complete failure to 



form a rotationally supported disk. The turbulent MHD 
model, on the other hand, shows good agreement with 
the Keplerian curve at least inside the radius of « 120 AU 
where the disk forms. Its rotation velocity profile is also 
very similar to the one of the MHD model with constant 
anomalous resistivity. Both models are able to reduce 
the magnetic braking effects by removing magnetic flux 
from the inner region and the resulting rotation curves 
of the formed disks are nearly Keplerian. In the resistive 
model, this is provided by the hyper-resistivity, while in 
the turbulent model is the turbulent reconnection that 
provides this diffusion. 

The bottom right panel of Figure 3 compares the pro- 
flies of the vertical component of the magnetic field, i?z, 
in the equator of the four models of Figure 2. While 
in the ideal MHD model, the intensity and gradient of 
the magnetic field in the central regions are very large 
due to the inward advection of magnetic flux by the col- 
lapsing material, in the anomalous resistive MHD model, 
the magnetic flux excess is completely removed from the 
central region resulting a smooth radial distribution of 
the fleld. In the turbulent model, the smaller intensity 
of the magnetic fleld in the inner region and smoother 
distribution along the radial direction compared to the 
ideal MHD case are clear evidences of the transport of 
magnetic flux to th e outskirts of the disk d ue to tur- 
bulent reconnection (|Santos-Lima et aI]|2010D . We note 
however that, due to the complex structure which is still 
evolving, the standard deviation from the average value 
is very large in the turbulent model, with a typical value 
of 100/iG (and even larger for radii smaller than 100 AU) 
which accounts for the turbulent component of the field. 

Finally, the bottom left panel of Figure 3 shows the 
mass of the formed disks in the four models, as a func- 
tion of the radius. In the hydrodynamical and the MHD 
resistive models {hydro and resistive) , the mass increases 
until R K, 250 AIJ and k, 350 AU, respectively, and both 
have similar masses. The masses in the ideal MHD and 
the turbulent MHD models {ideal MHD and turbulent) 
increase up to i? w 150 AU and k, 250 AU, respectively, 
and are smaller than those of the other models. Nonethe- 
less the turbulent MHD disk has a total mass three times 
larger than that of the ideal MHD model. 

4.1. Discussion 
4.2. Our approach and alternative ideas 

IShu et aO (|2006l ) (and references therein) mentioned 
the possibility that the ambipolar diffusion can be sub- 
stantially enhanced in circumstellar disks, but did not 
consider t his as a viable solution. The subsequent pa- 
per of Sh u et aP ([2007) refers to the anomalous resis- 
tivity and sketches the picture of magnetic loops be- 
ing reformed in the way of eventual removing magnetic 
flux. The latter process requires fast reconnection and 
we claim that in the presence of fast reconnection a more 
natural process associated with turbulence, i.e. magnetic 
reconnection diffusion can. solve the probl e m. 

Krasnopolsk v et all ([MfOi) and lLi et all (|2011[ ) showed 
by means of 2D simulations that an effective magnetic 
resistivity rj > 10^^ cm^ s~^ is needed for neutralizing 
the magnetic braking and enable the formation of a sta- 
ble, rotationally supported, 100 AU-scale disk around a 
protostar. The origin of this enhanced resistivity is com- 
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Fig. 2. — Face-on (top) and edge-on (bottom) density maps of the central slices of the collapsing disk models listed in Table 1 at a time 
t = 9 X IQil s (fs 0.03 Myr). The arrows in the top panels represent the velocity field direction an those in the bottom panels represent 
the magnetic field direction. From left to right rows it is depicted: (1) the pure hydrodynamic rotating system; (2) the ideal MHD model; 
(3) the MHD model with an anomalous resistivity lO"^ times larger than the Ohmic resistivity, i.e. iq = 1.2 X 10^" cm^ s~^; and (4) the 
turbulent MHD model with turbulence injected from t = until t=0.015 Myr. All the MHD models have an initial vertical magnetic field 
distribution with intensity = 35 /xG. Each image has a side of 1000 AU. 



pletely unclear and the value above is at least two to three 
orders of magnitude larger than the estimated ohmic dif- 
fusivity for these cores (e.g., [Krasnopolskv et al. 2010). 
On the other hand, these same authors found that am- 
bipolar diffusion, the mechanism often invoked to remove 
magnetic flux in star forming regions, is unabl e to pro- 
vide such required levels of diffusivity (see also iLi et al.l 

Eonl). 

In this work, we have explored a different mechanism to 
remove the magnetic flux excess from the central regions 
of a rotating magnetized collapsing core which is based 
on magnetic reconnection diffusion in a turbulent flow. 
Unlike the Ohmic resistivity enhancement, reconnection 
diffusion does not appeal to any hypothetical processes, 
but to the turbulence existing in the system and fast 
magnetic reconnection of turbulent magnetic flelds. One 
of the consequences of fast reconnection is that, unlike 
resistivity, it conserves magnetic fleld helicity. This may 
be important for constructing self-consistent models of 
disks. 

In order to compare our turbulent MHD model with 
other rotating disk formation models, we also performed 
3D simulations of a pure hydrodynamical, an ideal MHD 
and a resistive MHD model with a hyper-resistivity co- 
efficient ?7 ~ 10^° cm^ s~^ (Figure 2). The essential fea- 
tures produced in these three models are in agreement 



with the 2D models of iKrasnopolskv et al.l ()2010t ). i.e., 
the ideal MHD model is unable to produce a rotationally 
supported disk due to the magnetic flux excess that ac- 
cumulates in the central regions, while the MHD model 
with artificially enhanced resistivity produces a nearly- 
Keplerian disk with dimension, mass, and radial and ro- 
tational velocities similar to the pure hydrodynamical 
model. 

The rotating disk formed out of our turbulent MHD 
model exhibits rotation velocity and vertical magnetic 
fleld distributions along the radial direction which are 
similar to the resistive MHD model (Figure 3). These 
similarities indicate that the turbulent magnetic recon- 
nection is in fact acting to remove the magnetic flux ex- 
cess from the central regions, just like the ordinary en- 
hanced resistivity does in the resistive model. We note, 
however, that the disk formed out of the turbulent model 
is slightly smaller and less massive than the one produced 
in the hyper-resistive model. In our tests the later has a 
diameter ^ 250 AU, while the disk formed in the turbu- 
lent model has a diameter ~ 120 AU which is compatible 
with the observations. 

The effective resistivity associated to the MHD tur- 
bulence in the turbulent model is approximately given 
by -qturb ~ VturbLinj, whcre Vturb is the turbulent rms 
velocity, and Linj is the scale of injection of the turbu- 
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Fig. 3. — Radial profiles of the: (i) radial velocity vji (top left), (ii) rotational velocity (top right); (iii) inner disk mass (bottom left); 
and (iv) vertical magnetic field Bz, for the four models of Figure 2 at time t Ri 0.03 Myr). The velocities were averaged inside cylinders 
centered in the protostar v^ith height h = 400 AU and thickness dr = 20 AU. The magnetic field values were also averaged inside equatorial 
rings centered in the protostar. The standard deviation for the curves are not shown in order to make the visualization clearer, hut they 
have typical values of: 2 — 4 X 10* cm s~^ (for the radial velocity), 5 — 10 X 10* cm s^^ (for the rotational velocity), and 100 fiG (for the 
magnetic field). The vertical lines indicate the radius of the sink accretion zone. 



lence ^. We have adjusted the values of Linj and Vturb 
in the turbulent model employing turbulent dynamical 
times {Linj /Vturb) large enough to ensure that the cloud 
would not be destroyed by the turbulence before form- 
ing the disk. We tested several values of rjturt and the 
one employed in the model presented in Figure 2 is of 
the same order of the magnetic diffusivity of the model 
with enhanced resistivity, i.e., ijturb ^ VturbLinj ~ 10^" 
cm^ s~^. Smaller values were insufficient to produce rota- 
tionally supported disks. Nonetheless, further systematic 
parametric study shoul d be performed in th e futur e. 

As mentioned in §1, Mac hida etldl ()2010l . 12011] ) have 
also performed 3D MHD simulations of disk formation 
and obtained a rotationally supported disk solution when 
including only Ohmic resistivity (with a dependence on 
density and temperature obtained frorn the fi tting of the 
resistivities computed in lNakano et al.ll2002D . However, 
they had to evolve the system much longer, about four 
times longer than in our turbulent simulation, in order 
to obtain a rotationally supported disk of 100- AU scale. 

® We note however, that this value may be somewhat larger in 
the presence of the gravitational field (see L05, SXIO) 



In their simulations, a tiny rotationally supported disk 
forms in the beginning because the large Ohmic resistiv- 
ity that is present in the very high density inner regions 
is able to dissipate the magnetic fields there. Later, this 
disk grows to larger scales due to the depletion of the 
infalling envelope. Their initial conditions with a more 
massive gas core (which has a central density nearly ten 
times larger than in our models) probably helped the for- 
mation of the rotating massive disk (which is almost two 
orders of magnitude more massive than in our turbulent 
model). The comparison of our results with theirs in- 
dicate that even though at late stages Ohmic, or more 
possibly ambipolar diffusion, can become dominant in 
the high density cold gas, the turbulent diffusion in the 
early stages of accretion is able to form a light and large 
rotationally supported disk very quickly, in only a few 
10^ yr. 

Finally, we should remark that other mechanisms to re- 
move or reduce the effects of the magnetic braking in the 
inner regions of protostellar cores have been also inves- 
tigate d in the literature recently. iHennebelle &: Ciardil 
(|2009D verified that the magnetic braking efficiency may 
decrease significantly when the rotation axis of the core 
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is misaligned with the direction of the regular magnetic 
field. They claim that even for small angles of the order 
of 10 — 20° there are significant differences with respect 
to the aligne d case. Also, in a concomi tant work to the 
present one, iKrasnopolskv et al.l ()2011[ ) have examined 
the Hall effect on disk formation. They found that a 
Hall-induced magnetic torque can diffuse magnetic flux 
outward and generate a rotationally supported disk in 
the collapsing flow, even when the core is initially non- 
rotating, however the spun-up material remains too sub- 
Keplerian (Li et al. 2011). 

Of course, in the near future, these mechanisms must 
be tested along with the just proposed turbulent mag- 
netic reconnection and even with ambipolar diffusion, 
in order to assess the relative importance of each ef- 
fect on disk formation and evolution. Nonetheless, since 
MHD turbulence is expected to be present in these mag- 
netic cores (e.g., J 3allestcros-Parcdcsfc Mac Low 120021: 
MeUoli et all l2ci06nLeao et al. 2009; Sa ntos-Lima et al.l 
20101 and references therein), turbulent reconnection 
arises as a natural mechanism for removing magnetic flux 
excess and allowing the formation of these disks. 

4.3. Present work and SXIO 

In recent numerical study in SXIO, we showed that 
magnetic reconnection in a turbulent cloud can efficiently 
transport magnetic flux from the inner denser regions to 
the periphery of the cloud thus enabling the cloud to 
collapse to form a star. 

Here, also by means of fully 3D MHD simulations, we 
have investigated the same mechanism acting in a rotat- 
ing collapsing cloud core and shown that the magnetic 
flux excess of the inner regions of the system can be effec- 
tively removed allowing the formation of a rotationally 
sustained protostellar disk. 

Another empirical finding in SXIO is that the efficiency 
of the magnetic field expulsion via reconnection diffusiv- 
ity increases with the source gravitational field. This is 
a natural consequence of diffusion in the presence of the 
gravitational field which pulls one component (gas) and 
does not act on the other weightless component (mag- 
netic field). In terms of the problem in hand, this implies 
that more massive protostars can induce magnetic field 
segregation even for weaker level of turbulence. We plan 
to explore numerically this issue in a forthcoming work. 

4.4. Present result and bigger picture 

In this paper we showed that the concept of reconnec- 
tion diffusion (LOS) successfully works in the formation 
of protostellar disks. Together with our earlier testing 
of magnetic field removal through reconnection diffusion 
from collapsing clouds this paper supports a considerable 
change of the paradigm of star formation. Indeed, in the 
presence of reconnection diffusion, there is no necessity 
to appeal to ambipolar diffusion. The latter may still 
be important in low ionization, low turbulence environ- 

We note, however, that in the present study, we were in a highly 
magnetized disk regime, where the magneto-rotational instability 



ments, but, in any case, the domain of its applicability 
is seriously challenged. 

The application of the reconnection diffusion concept 
to protostellar disk formation and, in a more general 
framework, to accretion disks in general, is natural as 
the disks are expected to be turbulent, enabling our ap- 
peal to LV99 model of fast reconnection. An impor- 
tant accepted source of turbulence in accretion disks 
is the well known magneto-rotational instability (MRI) 
(Chandrasekhar 1960, Balbus & Hawley 1991) , but 
at earlier stages turbulence can be induced by the hy- 
drodynamical motions associated with the disk forma- 
tion. Turbulence is ubiquitous in astrophysical envi- 
ronments as it follows from theoretical considerations 
based on the high Reynolds numbers of astrophysi- 
cal flows and is strongly supported by studies of spec- 
tra of the interstellar electron density fluctuations (see 
[Armstrong et al.l 119951: [ Chepurnov fc LazarianI 120101 ) as 
well as of HI fl^azari anl 120091 for a review and refer- 
ences there in; 'Chcp urnov et al.ll2010[ ) and CO lines (see 
iPadoan et al. 2009). The application of the reconnection 
diffusion mechanism to already formed accretion disks 
will be investigated in detail elsewhere. It should be 
noted however that former studies of the injection of tur- 
bulence in accretion disks have shown that at this stage 
turbulenc e may be ineffective to mag netic flux diffusion 
outward (|Rothstein fc Lovelacell2008[ ). 

5. SUMMARY 

Appealing to the LV99 model of fast magnetic recon- 
nection and inspired by the successful demonstration of 
removal of magnetic field through reconnection diffusion 
from numerical models of molecular clounds in SXIO we 
have performed numerical simulations and demonstrated 
that: 

1 . The concept of reconnection diffusion is applicable 
to the formation of protostellar disks with radius ^100 
AU. The extension of this concept to accretion disks is 
foreseen. 

2. In the gravitational field, reconnection diffusion mit- 
igates magnetic breaking allowing the formation of pro- 
tostellar disks. 

3. The removal of magnetic field through reconnection 
diffusion is fast enough to explain observations without 
the necessity of appealing to enhanced fluid resistivity. 
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